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Abstract 



The Brueckner-Hartree-Fock formalism is applied to study spin polar- 
ized neutron matter properties. Results of the total energy per particle as a 
function of the spin polarization and density are presented for two modern 
realistic nucleon-nucleon interactions, Nijmegen II and Reid93. We find that 
the dependence of the energy on the spin polarization is practically parabolic 
in the full range of polarizations. The magnetic susceptibility of the system is 
computed. Our results show no indication of a ferromagnetic transition which 
becomes even more difficult as the density increases. 
PACS numbers: 26.60. +c, 21.60.Jz, 26.50.+X 



I. INTRODUCTION 

The study of the magnetic properties of dense matter is of considerable interest in con- 
nection with the physics of pulsars. These objects, since the suggestion of Gold IQ, are 
generally believed to be rapidly rotating neutron stars with strong surface magnetic fields 
of the order of 10^^ Gauss. Several authors have studied the possible existence of a phase 
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transition to a ferromagnetic state on pure neutron systems at densities corresponding to 
the theoretically stable neutron stars. Brownell and Callaway 0, and Rice considered a 
hard sphere gas model and showed that the ground state of the neutron gas becomes ferro- 
magnetic at fci? ~ 2.3 fm~^. Silverstein |Q and 0stgaard found that the inclusion of long 
range attraction significantly increased the ferromagnetic transition density {e.g. 0stgaard 
predicted the transition to occur at kp ~ 4.1 fm~^ using a simple central potential with 
hard core only for singlet spin states). Clark and Pearson and Saunier calculated 
the magnetic susceptibility for low densities {kp < 2 fm^^) using more realistic interactions. 
Pandharipande et al. using the Reid soft-core potential, performed a variational calcu- 
lation arriving to the conclusion that such a transition was not to be expected for kp < 5 
fm^^. Early calculations of the magnetic susceptibility within the Brueckner theory were 
performed by Backmann and Kallman employing the Reid soft-core potential, and re- 
sults from a Correlated Basis Function calculation were obtained by Jackson et al. [jlO[ with 



the Reid Vq interaction. A different point of view was followed by Vidaurre et al. |TT|, who 
employed neutron-neutron effective interactions of Skyrme type, finding the ferromagnetic 
transition at kp ~ 1.73-1.97 fm~^. 

In connection with the problem of the neutrino diffusion in dense matter, Fantoni et al. 



I4 have recently employed a new quantum simulation technique (the so-called Auxiliary 
Field Diffusion Monte Carlo method (AFDMC)) using realistic interactions (based upon the 
Argonne vig two-body potential [|1^ plus Urbana IX three-body potential [0), and have 



found that the magnetic susceptibility of neutron matter shows a strong reduction of about 
a factor 3 with respect its Fermi gas value. They pointed out that such a reduction may 
have strong effects on the mean free path of a neutrino in dense matter and, therefore, it 
should be taken into account in the studies of supernovae and proto-neutron stars. 

In this work we employ the Brueckner-Hartree-Fock (BHF) approximation, using the 
realistic Nijmegen II and Reid93 |[T^ nucleon-nucleon interacions, to study spin polarized 
neutron matter properties such as the total energy per particle and the magnetic suscepti- 
bility. We employ the so-called continuous prescription when solving the Bethe-Goldstone 



equation. As shown by Song et al. [|T5| , the effects from three-body clusters are dimished in 
this prescription. We also explore in this work the dependence of the total energy per parti- 
cle on the spin polarization, finding that up to a very good approximation this dependence 
is parabohc. 

The paper is organized in the following way. In Section |I| the theoretical background 
of our calculation is briefly reviewed. The construction of the neutron-neutron G-matrices 
and the calculation of the total energy per particle are shown in Section |II A| , whereas the 
magnetic susceptibility is determined in Section |IIB| . Our results are presented in Section 
|T|. Finally, a short summary and the main conclusions are given in Section 



II. THEORETICAL BACKGROUND 

In this section we briefly show how to evaluate, in the BHF approximation, the total 
energy per particle and the magnetic susceptibility of a system of neutrons in which we 
assume that the density of particles with spin up, p\ is different from that with spin down, 

A. Energy per particle 



Our calculation of the total energy per particle starts with the construction of G-matrices, 
which describe in an effective way the interaction between two neutrons, for each one of the 
spin combinations (It) Ti? it or ||), in the presence of a surrounding medium. They can 
be obtained by solving the following integral Bethe-Goldstone equations 
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In the above expressions the first (last) two subindices indicate the spin projections of the 
two neutrons in the initial (final) state, V is the bare nucleon-nucleon interaction, Q is the 
Pauli operator which allows only intermediate states compatible with the Pauli principle, and 
u is the starting energy defined as the sum of single-particle energies, et{i)' '^^ interacting 
neutrons. Note that G^i^^i and Gji^^i are obtained from a coupled channel equation due to 
the mixing induced by the interaction. One can equivalently solve the Bethe-Goldstone 
equation in the spin-coupled basis, where the interaction is diagonal, although in that case 
the Pauli operator is non-diagonal. However, the conventional angle-average of the Pauli 
operator makes it diagonal, thus reducing the problem to an uncoupled one in each total 
spin channel. 

The single-particle energy of a neutron with momentum k and spin projection a =t(i) 
is given by 

where the single-particle potential Ua{k) represents the average field felt by the neutron due 
to its interaction with the other neutrons of the system. In the BHF approximation it is 
given by 

U^ik) = Re J2 E {kk'\G^^>,^^itu = e, + e,.)\kk')A , (3) 
o"'=T,J. fc'<fcj' 

where a sum over the two Fermi seas of spin up and down, characterized by k\ = {Gir'^p^y^^ 
and kp = (Gvr^pJ-)^/^ respectively, is performed and the matrix elements are properly anti- 
symmetrized. 

Once a self-consistent solution of Eqs. (|l]) and (^ is obtained, the total energy per 
particle is easily calculated 

This quantity is a function of and or, equivalently, of the total density p = + and 
the spin polarization A, defined as 
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A = ^. (5) 

Note that the value A = corresponds to non-polarized or paramagnetic (p^ = p^) neutron 
matter, whereas A = ±1 means that the system is totally polarized, i.e., all the spins are 
aligned in the same direction. 

B. Magnetic susceptibility 

The magnetic susceptibility of a system characterizes the response of this system to a 
magnetic field and gives a measure of the energy required to produce a net spin alignment 
in the direction of the field. It is defined as 

(dM\ 

where M. is the magnetization of the system per unit volume given by 

M = /x(pT - p^) = ppA , (7) 

with the magnetic moment of a neutron, and Ti is the magnetic field which can be obtained 
from 

( d(E/N) \ 1 ( d{E/N) \ 

Using Eqs. (|^) and (H), the magnetic susceptibility can be written as 

y = (9) 

where the second derivative can be taken at A = if the field Ti. is assumed to be small. 

It is customary to study the magnetic susceptibility in terms of the ratio x/Xf, where 
Xf is the magnetic susceptibility of a free Fermi gas, usually known as Pauli susceptibility. 
It can be straightforwardly obtained from Eq. (|^) and the total energy per particle of the 
free Fermi gas 
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Xf = TT^kp . 10) 

where the Fermi momentum kp = {Stt^pY^^ is related to kp and kp through the relations 

k^p = kp{i + Ay/' ^^^^ 



III. RESULTS 

The total energy per particle for totally polarized (solid lines) and non-polarized (dashed 
lines) neutron matter is shown in Fig. |I| as a function of the density. Results for the Nijmegen 
II interaction are plotted on the left panel, whereas those corresponding to the Reid93 
interaction are shown on the right panel. As can be seen from the figure, for both interaction 
models, the energy of totally polarized neutron matter is always more repulsive than non- 
polarized neutron matter in all the density range explored. This additional repulsion can 
be understood, firstly, in terms of the kinetic energy contribution, which is larger in the 
totally polarized case than in the non-polarized one; and secondly, in terms of the potential 
energy contribution because, due to symmetry arguments, all partial waves with even orbital 
angular momentum L (some of them attractive, as the important ^Sq) are excluded in totally 
polarized neutron matter. In order to illustrate this, we have plotted in Fig. ^ the separate 
kinetic (left panel) and potential (right panel) energy contributions for the Nijmegen II 
interaction model (similar results are obtained for the Reid93 one, but they are not included 
in order to make the discussion more clear). An interesting conclusion which can be inferred 
from these results is that a phase transition to a ferromagnetic state is not to be expected 
from our calculation. If such a transition would exist a crossing of the energies of the totally 
polarized and the non-polarized cases would be observed at some density, indicating that 
the ground state of the system would be ferromagnetic from that density on. As can be seen 
in Fig. |I|, there is no sign of such a crossing and, on the contrary, it becomes less favourably 
as the density increases. 
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We have shown results for totally polarized and non-polarized neutron matter. Let us 
consider now an intermediate situation in which not all the spins, but a part of them, are 
aligned in a given direction, and let us examine the dependence of the total energy per 
particle in the spin polarization A. This dependence is shown in Fig. ^ for five different 
densities {po/2, pQ,2pQ,5po and 7po, being po = 0.17 fm~^ the saturation density of nuclear 
matter). As in Fig. |l], results for the Nijmegen (Reid93) are shown on the left (right) 
panel. Circles, squares, diamonds and triangles correspond to our BHF results, whereas 
solid lines correspond to the parabolic approximation discussed below. As can be seen from 
this figure, and as it was expected, E/N is symmetric in A. It can also be seen in this figure 
that E/N shows a minimum at A = for all the densities considered, being this again an 
indication that the ground state of neutron matter is paramagnetic. Another interesting 
thing is to note that this dependence is up to a very good approximation parabolic, being 
this parabolic character only slightly lost at very large densities. Therefore, in the same 
spirit as it is done in nuclear matter to determine the symmetry energy, one can try to 
characterize this dependence in the following simple analytic form 

^(p,A) = |(p,0) + a(p)A% (12) 

where, assuming the quadratic dependence to be valid up to | A |= 1 as our results indicate, 
the value of a(p) can be easily obtained for each density as the difference between the total 
energy per particle of totally polarized and non-polarized neutron matter 

a(p) = ^(p,±l)-^(p,0) . (13) 

The magnetic susceptibility can be evaluated in a very simple way if the parabolic de- 
pendence of Eq. (|1^) is assumed, giving 

^(-' - ^) ■ 

In Fig. 1^ the ratio x/Xf is shown as a function of the density. The solid line shows the 
result for the Nijmegen II interaction, whereas the dashed line corresponds to the one ob- 
tained with Reid93. Starting from 1, the ratio decreases rapidly for small densities and more 



slowly as density increases. It can be inferred again from this figure that a ferromagnetic 
phase transition, which would be signaled by an infinite discontinuity giving rise to a change 
of sign in xIxf-, is not seen and not expected at larger densities either. 

Finally, our results for x/Xf are compared in Table | with those of the recent calculation 
performed by Fantoni et al. [0, shown in columns labelled AU6' and AU8'. As can be 
seen from the table, there is a very good agreement between these results and ours. For 
completeness, we show in parentheses the results obtained when the standard discontinuous 
prescription is used in solving the Bethe-Goldstone equation. In both prescriptions, the 
results are very similar, which is not surprising due to the fact that x is obtained from an 
energy difference (see Eqs. (|13D and(0)) which partly cancels the possible discrepancies. 
Only for densities larger than 2po, the discontinuous prescription results differ more than 
10% from the continuous ones. 



IV. SUMMARY AND CONCLUSIONS 

Employing realistic modern nucleon-nucleon interactions (Nijmegen II and Reid93) we 
have performed a Brueckner-Hartree-Fock calculation of spin polarized neutron matter 
properties. We have studied the total energy per particle of neutron matter as a func- 
tion of the density and the spin polarization A. We have found that in the range of densities 
explored (up to 7po) totally polarized matter is always more repulsive than non-polarized 
matter, being this an indication that a phase transition of the system to a ferromagnetic 
state is not expected. 

We have seen that the total energy per particle is not only symmetric on the spin polar- 
ization A, as it was expected, but also parabolic in a very good approximation up to | A |= 1 
even at high densities. This finding supports the calculation of the magnetic susceptibility of 
neutron matter by using only the energies of the spin symmetric and fully polarized systems. 

Finally, we have calculated the magnetic susceptibility of the system as a function of the 
density, finding a very good agreement with a recent Monte Carlo calculation |]T2 . 
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TABLES 

TABLE I. Magnetic susceptibility ratio x/xf- Our BHF results, labelled Nijmegen II and 
Rcid93, are compared with the AFDMC results of Fantoni et al. [12], labelled AU6' and AU8'. 
BHF results obtained with the standard discontinuous prescription are given in parentheses. 



p/po 


Nijmegen II 


Reid93 


AU6' 


AU8' 


0.75 


0.39(0.41) 


0.39(0.41) 


0.40 




1.25 


0.38(0.39) 


0.37(0.39) 


0.37 


0.39 


2.0 


0.34(0.37) 


0.34(0.38) 


0.33 


0.35 


2.5 


0.32(0.36) 


0.33(0.37) 


0.30 
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FIG. 1. Total energy per particle as a function of the density for totally polarized (solid lines) 
and non-polarized (dashed lines) neutron matter. The left panel shows results for the Nijmegen II 
nucleon-nucleon interaction, whereas results on the right panel correspond to the Reid93 interac- 
tion. 
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FIG. 2. Kinetic (left panel) and potential (right panel) energy contributions to the total energy 
per particle as a function of the density for totally polarized (solid lines) and non-polarized (dashed 
lines) neutron matter. Results are shown for the Nijmegen II interaction. 
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FIG. 3. Total energy per particle as a function of the spin polarization A for different densities. 
Results for the Nijmegen II (Reid93) interaction are shown in the left (right) panel. Circles, 
squares, diamonds and triangles show our BHF results, whereas solid lines correspond to the 
parabolic approximation defined in Eq. (p^. 
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FIG. 4. Ratio x/xf as a function of the density. The sohd hne shows the result for the Nijmegen 
II interaction, while the dashed line corresponds to the one obtained with Reid93. 
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